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Gravity  Vector  Determination  from  Inertial  and  Auxiliary  Data  and 
Potential  Utilization  of  Generated  Vector  Component  Information 


H.  BAUSSUS  VON  LUETZOW,  DR. 

U.S.  Army  Engineer  Topographic  Laboratories 
Port  Belvoir,  Virginia  22060 


1.  INTRODUCTION.  The  present  Rapid  Geodetic  Survey  System  (RGSS), 
developed  by  Litton  Systems,  Inc.  for  the  U.S.  Army  Engineer  Topographic 
Laboratories  (ETL),  has  the  capability  of  determining  horizontal  and 
vertical  positions,  deflections  of  the  vertical,  and  gravity  anomalies 
vrith  average  rms  errors  of  1  m,  0.3  m,  2  arcsec,  and  2‘mgal,  respective¬ 
ly,  for  50  tan  runs  under  utilization  of  post-mission  adjustments.  It 
operates  as  a  quasi  local-level  system.  It  thus  does  not  require  alti¬ 
tude  damping,  permits  Kalman  stochastic  error  control  without  great 
complexity  under  consideration  of  observed  velocity  errors  at  vehicle 
stops,  and  employs  effective  post-mission  adjustments  with  the  aid  of 
terminal  position  and  gravity  vector  information  only.  Present  critical 
hardware  consists  of  an  A-1000  vertical  accelerometer,  two  A-200 
horizontal  accelerometers,  and  two  0-300  gyroscopes.  Although  data 
utilization  from  repeated  runs  has  resulted  in  average  deflection  rms 
errors  of  about  1  arcsec,  the  requirement  of  a  maximum  rms  error  of  0.5 
arcsec  necessitates  the  incorporation  of  higher  performance  gyroscopes 
and  accelerometers.  This  would  simultaneously  achieve  conventional 
surveying  accuracy  (10”°).  following  the  approximate  elimination  of  the 
effects  of  constant  gyro  biases  on -positions  under  utilization  of 
accurate  initial  and  terminal  coordinates,  the  remaining  dominant  error 
source  is  gyro  correlated  random  noise.  In  this  respect,  the  pertinent 
parameters  for  the  G-300  instruments  are  a  standard  deviation  of  0.002° 
hr”  and  a  correlation  time  of  3  hours.  Another  error  source  to  be 
considered  is  the  accelerometer  scale  factor,  assumed  to  be  constant  for 
a  test  run.  The  standard  scale  factor  error  for  the  A-200  accelerometer 
is  0.01*?.  Correlated  accelerometer  noise  is  characterized  by  a  standard 
deviation  of  10  mgal  and  a  correlation  time  of  40  minutes.  Both  corre¬ 
lated  gyro  random  drift  and  correlated  accelerometer  measurement  errors 
affect  the  accuracy  of  deflections  of  the  vertical  estimated  from 
inertial  data  and  initial  and  terminal  deflection  compdnents,  possibly 
augmented  by  corresponding  geodetic  azimuths.  In  order  to  achieve 
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deflection  accuracies  smaller  than  or  equal  to  0.5  arcsec  rms,  it  is, 
therefore,  necessary  to  install  higher  performance  gyros  and  accelero¬ 
meters.  A  significant  reduction  of  the  two  autocorrelation  parameters 
is  also  a  necessary  prerequisite  for  the  optimal  statistical  estimation 
of  deflections  of  the  vertical.  The  Litton  G-1200  gyro  is  expected  to 
be  compatible  with  the  stringent  deflection  accuracy  requirement.  Its 
random  drift  error  is  about  0.001°  hr"1  and  its  autocorrelation  parame¬ 
ter  is  short.  For  details,  reference  is  made  to  Litton  [l]  and  Huddle 
[2].  The  A-1000  accelerometer  addressed  by  Litton  [3*4-]  has  a  standard 
scale  factor  error  of  approximately  0.005!?.  Correlated  accelerometer 
noise  is  2  mgal  rms,  equivalent  to  about  0.35  arcsec. 

2.  DETERMINATION  OF  VERTICAL  GRAVITY  COMPONENT.  In  a  local  Cartesian 
coordinate  system  tangent  to  a  normal  gravity  surface,  with  x  as  west- 
east  coordinate,  y  as  south-north  coordinate  and  z  as  vertical  coordi¬ 
nate,  positive  upward,  it  is 

g  ■  (gn,  -  g€,  -g)  (1) 

In  eq.  (l),  g  is  the  gravity  vector  at  the  origin,  -g  is  its  vertical 
component,  and  gn  and  -g£  are  its  horizontal  components,  n  and  £  being 
the  prime  and  meridian  deflections,  respectively.  For  computational 
and  estimation  purposes  it  is  advantageous  to  use  the  gravity  anomaly  Ag 
■  g  -Y  with  y  as  normal  gravity,  to  be  found  by  consideration  of  the 
elevation  above  the  normal  reference  surface. 

The  RGSS  error  differential  equation  for  the  vertical  channel, 
including  first  and  second  order  time  derivatives,  is 

..9  1 


2o,a  z  -  (pN 


2  V  *  "  PE  y 


Hg  ♦  a 


+  A  5  -  A  6 

'  x  zx  y  zy 


where  the  symbol  6  of  position  errors  5x,  <5y,  6  z  is  omitted  for  simplic¬ 
ity.  In  eq.  (2),  ■  g  R  with  R  as  the  earth's  mean  radius  is  the 

squared  Schuler  frequency,  a  V  R"  with  V  as  the  system’s  east 

velocity  is  the  north  angular  rate,  is  the  north  earth  rate,  Pg  * 

-  V^R"  with  Vy  as  the  system's  north  velocity  is  the  east  angular  rate, 
az  is  the  vertical  accelerometer  error,  Ax  and  Ay  are  east  and  north 
accelerometer  outputs,  and  <5XZ  and  5yZ  are  z-axi3  misalignment  errprs  in 
the  x,z  and  y,z  planes.  At  the  system's  stop,  the  term  involving  x  may 
be  estimated  since  x>  is  observed,  and  the  last  two. terms  are  omitted. 

At  terrestrial  vehicle  speeds,  the  term  involving  y  becomes 
insignificant.  1 


See  Litton 


3-132. 
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Because  of  initial  calibration  z(tQ)  *  zQ  -  0  and  zQ«  iQ-  f0  -  Ajg 

••  •• 

*  A  *  0,  the  term  AgQ  +  az  has  to  be  subtracted  from  z(ty)  ■  Zy.  At 
a  survey  vehicle  stop,  therefore,  under  omission  of  pw±  and  p*#  in  eq. 
(2), 


Agy  -  Ago  ♦  zv  ♦  2nNxv  -  2aig  zy  -  (a 


z  -  az 
V  0 


(3) 


where  the  error  terms 


Zv  and  a^ 


a  are  weakly  correlated.  If  z 

o 


and  az  are  considered  as  random  variables  with  zero  means  and  estimated 
variances  and  covariances,  respectively,  the  gravity  anomaly  error 
variance  is 


var  Sg^ 


4 

4-oj  var  z 
s  v 


var  (a 

z 

v 


(4) 


Since  zy  has  an  average  rms  error  of  0.3  m  due  to  Kalman  filter  estima¬ 
tion  by  means  of  observed  zy  at  vehicle  stops  and  a  post-mission  adjust¬ 
ment  under  utilization  of  a  terminal  zp,  the  rms  error  induced  by  Zy  is 
only  0,1  mgal.  High  accuracy  determinations  of  gravity  anomalies  or  of 
the  vertical  component  of  gravity  by  means  of  a  local-level  inertial 
system  thus  require  a  high  performance  vertical  accelerometer,  high 
accuracy  height  determinations  and,  possibly,  utilization  of  i-observa- 
tions  at  vehicle  stops.1 

The  availability  of  Age  and  zQ  at  the  end  of  the  survey  permits  the 
determination  of  az  -  az  *  fp  in  eq.  (3).  Empirical  corrections  for 

a,  -  a_  in  the  form 


C  f 


v  e 


(5) 


prevent  error  growths  to  about  2  var  az«  Under  restrictive  assumptions, 
with  the  omission  of  the  subscript  z. 


c  3  cov  (a 

V  V 


I 


(6) 


3.  SINGLE  HORIZONTAL  CHANNEL  DETERMINATION  OF  DEFLECTIONS  OF  THE 
VERTICAL.  While  the  determination  of  gravity  anomalies  ‘may  be  less 
rapidly  accomplished  with  high  accuracy  by  means  of  gravimeter  and 
height  measurements,  astrogeodetic  and  gravimetric  determinations  of 
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deflections  of  the  vertical  are  time-consuming.  According  to  Moritz 
[6],  astrogeodetic  and  gravimetric  rms  errors  are  approximately  0.2 
arcsec  and  1.5  arcsec,  respectively.  Rapid  inertial-astronomic 
determinations  of  n  and  ?  with  rms  errors  not  exceeding  0.35  arcsec 
would  thus  be  highly  economical  and  satisfactory  for  many  purposes. 

Both  the  single  and  coupled  horizontal  channel  determinations  of  n  and  £ 
require  the  utilization  of  the  simplified  error  differential  equations2! 


Vz  *  +  4 


(7) 


7  -  ~Ss*z  «■  g^s  -  g£  + 


(9) 


R  ^rafr.x  ♦  R  1  (fljj  ♦  sec  <j>)y 


*Ve  40  (9) 


r”1  x 


_d>  -c.  +6  (10) 

Z  i 


_  — 1  . 
-R  y 


*Vz 


n  (ii) 


Symbols  used  in  the  foregoing  equations  including  total  time  derivatives 
of  first  and  second  order  are,  except  for  those  defined  in  section  2, 


$z  azimuth  platform  attitude  error 

<tfl  platform  tilt  error  about  north  axis 

platform  tilt  error  about  east  axis 
east  acceleration  of  survey  vehicle 
S.J  north  acceleration  of  survey  vehicle 

a^  correlated  east  accelerometer  error 

correlated  north  accelerometer  error 
$  geographic  latitude 

(i»2  '  3  *  R-  tn  $  Vx  vertical  spatial  rate 

djj  *  *  p.j  north  spatial  rate 

nig  *  -H"  Vy  east  spatial  rate 

a  azimuth  axis  angular  drift  rate  error 

B  north  axis  angular  drift  rate  error 

Y  east  axis  angular  drift  rate  error 

For  land  vehicles,  eqs.  (9)  -  (ll)  may  be  simplified  by  omission  of  p^, 


2.  A  more  complex  system  may  be  required.  See  Litton- (5],  pp.  3-252  - 
3-254. 


4 


BAUSSUS  VON  LUETZOW 


ug,  and  p2  *  R“  tn$.Vx,  and  by  use  of  constant  accelerations  Sjj  and  Sg 
which  should  be  approximately  achieved.  Then,  u>z  *  Qsiafr  and  *ftcos$ 
where  fl  denotes  the  earth’s  inertial  angular  velocity.  The  initial 
conditions  at  t0»  0  are,  un$er  consideration  of  plumbline  leveling,  wz(0)  * 
0,  $5(0)  -  n0,  #E(0)  -  50,  x(o)  -  y(o)  -  x(o)  -  y(o)  -  0. 

In  the  general  RGSS  mode  of  operation,  position  and  gravity  vector 
information  is  generated.  The  Kalman  filter  mechanization  is  simplified 
by  setting  ug  *  »  0  in  eqs.  (9)  -  (ll)  which  results  in  the 

uncoupling  of  the  horizontal  channel  differential  equations.  This  is 
not  critical  in  the  context  of  present  instrumentation  for  travel 
periods  of  2  hours.  As  an  example,  eqs.  (3)  and  (10)  may  then  be  formu¬ 
lated  at  the  end  of  the  first  travel  period,  whence  t  =  t1  , 

7,  58  «<t>1  -  -  a>}  (12) 

1  ‘ov. 

1 

-  -H  7  +  /0  Ydt  +  '  (13) 

In  eq.  (12),  accelerometer  calibration  at  t  ®  0  ha3  been  taken  into 
account.  The  quantity. -R“  y  in  eq.  (13)  is  well  estimated  by  the  ob¬ 
served  velocity  error  y«  After  implementation  of  the  corresponding  tilt 
correction,  under  neglect  of  insignificant  estimation  errors  and  consid¬ 
eration  of  the  bias  error  <£g(c^), 


y(tj  -  3  -  €  (Ci  -  CQ)  +  g  I1  (y  -  y)  dt 

o 


+  a:).  -  aM 


(u) 


y(tn)  3  yn  3  gfB  -  g  (5n  -0 


O  ♦  g  fn  (y  -  y)  dt 


+  3n  -  aN 
n  o 


(15) 


The  propagation  of  constant  gyro  bi^s. errors  may  be  approximated  by 
a  solution  of  eqs.  (9)  -  ( 1  1 )  with  y  =  x=  y  =  0  and  by  treating 

and  ojg  as  constants  in  accordance  with  Huddle  [§].  The  result  is 


K  3  *  (r 


+  stn  /2  - 

z  N 


(16) 


t  (-  YtQ/2  ♦  8 
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) 


(17) 


3  t  (ytn^/2  ♦  a  )  (18) 

Linder  consideration  of  terminal  deflection  and  azimuth  errors  at  t  ■  t. 
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it  is  possible  to  estimate  5,  S ,  ?  from 


h 
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N 

nn 
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nn 

(19) 
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■  ml 

estimated 

A 

n 

observed 
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Equations  (15)  and  (16)  show  that  the  computation  of  fg  in  eq.  (14)  is 
to  the  first  order  associated  with  an  error  * 

t.  t 

"  — ~r~  /0  ^  dt  +  ^  (2°) 


The  final  solution  for  the  meridan  deflection  of  the  vertical  from  eqs. 
(14)  and  (20)  is 

7, 


5,  -  tin  *  -  --) 

l  o  Ei  g 


aN.  "  aN  t.  ®S  “  aN 

i  o  x  n  o 
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n  -  '  (21) 

♦  /q1  -  7  )dt  -  £  n  (y  -  Y  )dt 


Equation  (21)  may  be  supplemented  by 


t.  t . 

5C.  •  «?  +  -mi-  (6C_  -  Hj  -  (*„  -  -1  4~  ) 


n 


e 


o .  t  o 
i  n  n 


(22) 


to  account  for  astrogeodetic  deflection  errors  and  accelerometer  bias 
errors.  For  a  straight  traverse,  the  last  term  in  eq.  (22)  tends  to 
cancel  out. 


The  rms  deflection  error  Ojr(t^)  can  be  qoraputed  by  covariance  analysis 
involving  the  terras  without  parentheses  in  eq.  (21).  Under  inclusion  of 
the  first  two  terms  of  eq.  (22), 

t .  2  t .  2 

var  C.  =  var  a.  *  var  y .  +  r  1  -  )  var  5  +  (■—“)  var  ?  ,  (25) 

l  l  i't  o  t  n 

n  n 

where  var  a^  is  the  accelerometer-induced  variance  and  vary^  designates 
the  gyro-induced  variance.  Under  consideration  of  present  ROSS  param¬ 
eters,  identified  ih  section  1,  the  rms  deflection  error  a.  is  approxi¬ 
mately  presented  in  Figure  1  for  a  vehicle  speed,  of  50  tan/hi 
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? 


Distance  Traveled  (km) 

Figure  1 

Approximate  RMS  Defleccion  Error  as  a  Funccion  of  Travel  Distance 

The  variation  of  normal  gravity  g  with  altitude  is  not  critical  in 
terrestrial  applications.  The  scale  factor-generated  error  does  not 
vanish  with  respect  to  L-shaped  traverses.  It  is,  however,  relatively 
siuoll  for  travel  times  not  exceeding  2  hours.  Heading  sensitivity 
induced  by  significant  azimuth  changes  appears  to  be  somewhat  critical 
and  may  require  empirical  corrections.  Schwarz  [7]  discussed  the  impact 
of  heading  sensitivity  on  positioning  and  concluded  chat  error 
reductions  may  be  achieved  by  a  modified  smoothing  procedure.  The 
problem  is  more  intricate  in  the  case  of  deflection  determinations.  The 
use  of  static  accelerometer  measurements  in  conjunction  with  initial 
and  terminal  deflections  together  with  a  simplified  Kalman  filter  was 
first  discussed  by  Huddle  [8].  An  advanced  method  for  single  channel 
vertical  defleccion  determination  was  presented  by  Saussus,  von  Luetzow 
[9].  It  utilizes  all  observations  simultaneously  and  appears  in  a 
modified  form  in  section  4. 

4.  COUPLED  HORIZONTAL  CHANNEL  OPTIMAL  DETERMINATION  OF  DEFLECTIONS  OF 
THE  VERTICAL  IN  SEMI-FLAT  TERRAIN.  Higher  accuracy  requirements 
necessitating  installation  of  gryoscopes  and  accelerometers  with  small 
error  variances  and  short  correlation  times  as  performance  character¬ 
istics  require  che  integration  of  the  system  of  differential  equations 
for  besc  possible  estimation  of  n  and  ?  and  for  error  assessment  in  che 
case  of  single  channel  solutions.  For  this  purpose,  eqs.  (7)  and  (8) 
are  wricten  as 

X  "  SN*z  "  3,N  +  3(711  "  V  +  aE  “  aE  (24) 

o 

y  *  -S^  +  g*E  -  8(5  -  5q)  +  aN  -  aN  -  <  (25) 

o 

In  this  wav,  initial  accelerometer  calibrations  have  been  taken  inco 
account  and  the  initial  conditions  are 
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$s(o)  -  (o)  -  $z(o)  *  x(o)  *  7(0)  m  x(o)  ■  7(0)  -  0  (26) 


The  accelerations  Sjj  and  Sg  can  probabl7  be  neglected.  If  the  37stem  is 
treated  as  one  with  constant  coefficients  as  a  good  approximation,  a 
closed  solution  as  a  function  of  time  is  possible.  Because  of  intermit¬ 
tent  Kalman  filter  corrections  and  the  need  for  numerical  weight  factors 
j,t  appears  to  be  advantageous  to  attempt  a  numerical  solution  for  x  and 
7  under  utilisation  of  terminal  deflection  and  azimuth  data.  7or  eco¬ 
nomic  reasons  it  is  assumed  that  the  surve7  vehicle  travels  approxi¬ 
mate^  at  a  constant  speed  when  in  motion  and  stops  ever7  3  minutes  for 
1  minute.  The  speed  should  not  exceed  10  msec-  in  order  to  restrict  the 
length  of  travel  intervals. 


Solutions  for  $g,  z * 
the  integration  schemes 


x, 


x,  7  are  obtained  in  accordance  with 


F  -*■  F  At 
0  o 


.  (27) 


2F 


v-1 


-  F  + 
v-2 


(At) 


"  v  — 1 


(28) 


with  At  *  30  sec  and  possibl7  60  sec.  The  solution  structure  at  the  end 

of  the  first  stop  interval,  indicated  b7  the  subscript  s*1 ,  is 

v  v  v  v  v 

Tea  +TdB  ♦  |  e  T  •  (29) 

A  v  v  A  A  cv'v  v 


■\>1 


T  a  v  +  )  b  w  + 


v 

\l  V 


y  b  w 

^  v  v 


1 


In  eq.  (29),  v  *  g  (n  -  nQ)  ♦  ag 


V 


-g(5  5q) 


a«  -  a- 


H. 


Under  consideration  of  Kalman  filter  tilt  corrections,  assumed  for 
simplicity  to  eliminate  the  integrated  first  two  terras  in  eq.  (9)  and 
the  first  term  in  both  eq.  (10)  and  ( 1 1 7 » 


v 

l  A  v 
^  v  v 


V 


♦  T. 


(30) 


where  T^  and  subsequent  Tg  represents  an  aggregate  of  tilt-induced 
randoip  erpors.  During  the  stop  interval,  eqs.  (9)-(ll)  are  integrated 
with  x  *  7  -  7  -  0  with  a  resultant  effect  on  eqs.  (7)  and  (8).  The 
integration  is  then  continued  and  7ields  under  utilization  of  average 
biases  the  solution  structure 
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,(2)  .  ?• 


T3  Av+T  Bw+Ca+DB+Ey+R  +  T 
L  v  v  4  v  u  v  v  v  v  3  3 


where  R^  designates  residual,  random-type  terms  including  a  -  av ,  ff  - 

8V,  and  f  -  yv> 

At  the  termination  of  the  survey,  when  nn,  £n,  and  An  are  avail¬ 
able,  it  is  possible  to  solve  for  the  gyro  biases  in  the  form 


51  I  F 


3  3  Eg  +  e2  +  ^2 

.yj  tay  U^J  Uy 

where  the  P's  are  computable,  the  e's  are  an  aggregate  of  v  and  w-er- 
rors,  and  the  5's  are  aggregates  of  gyro  random  and  tilt- errors. 
Substitution  of  a,  S,  Y  from  eq.  (32)  in  eq.  (30  and  separation  of  v 
and  w  into  deflection  a$d  accelerometer  errors  leads,  tinder  restriction 
to  solutions  for  x  and  y,  to  the  final  results 

n  n  -1”(2)  -(2) 
n+y  k  n  +  £  1  £  *  n  +  g  xv  +  +  r  +  r  +  r  +■  r. 

'3  u  v  L  v  v  o  s  N  gs  as  ts  as 


n  n  _1  "f2")  -(2  ) 

g  *  y  mn  +  y  nC  *  £  -  g  y  +  <t„  '+  p  +  p  +■  p,  +  p,  (34-) 

*3  L  '  v  V  “v’v  ’0  s  's  yEs  gs  as  ts  ds  w 

u*0  v=*0 

The  last  4  terms  on  the  right  aide  of  both  eqs.  (33)  and  (34)  are  aggre¬ 
gates  of  random  errors  associated  with  gyro3,  accelerometers,  tilt 
corrections,  and  initial  and  terminal  deflection  errors.  Equations  (33) 
and  (34)  are  reformulated  as 
n  n 

n+y  kn+y  1  £  -  S  =  M  -  H  (35) 

s^-vv^-^v  v  n  n  n 
u=0  \>=0  333 


n  .  n 

£  ♦  y  mn  +y’n£ 

3  L  „  V  V  L  *  V  V 

v=0  vfO 


h  =  V  "  V 

333 


where  S,  M,  N  denote  signal,  measureable  message',  and  no'n-measurable 
noise,  respectively. 


fJnder  utilization  of  deflection  of  the  vertical  covariance  functions. 
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the  prime  deflection  of  the  vertical  can  be  optimally  estimated  in  the 
form 


ie 


]  +  B 


ie 


(37) 


where  i-s,  Aie  and  Bie  are  matrices  of  regression  coefficients,  and  the 
terms  in  brackets  are  message  matrices.  With  lc  -  0,1,...., n,  A^e  and 
8^q  can  be  determined  from  the  equations 


rsn 


e1-  n, 

L 


He 


+  B. 
le 


■V 


(38) 


A.  i_Sr  S„ 
18  5i  \ 


♦  N, 


iC 


*  V-3*. 


+  H, 


5 .  S5.  1 

1  lc 


(39) 


In  eqs.  (38)  and  (39),  where  the  bar  symbol  stands  for  covariance,  the 
noise  covariances  need  only  be  computed  once. 

Simplified  solutions,  particularly  in  the  case  of  approximately 
straight  traverses,  are  possible.  The  coupling  of  eqs.  (33)  and  (34) 
should,  however,  always  be  considered.  If  averaged  message-type  data 
from  repeated  surveys  are  employed,  the  instrument-generated  noise 
covariances  in  eqs.  (38)  and  (39)  are  to  be  reduced. 

5.  INTERPOLATION  OF  GRAVITY  ANOMALIES  AND  DEFLECTIONS  OF  THE  VERTICAL 
IN  MOUNTAINOUS  TERRAIN.*^  It  is  well  known  that  an  accurate  analytical 
representation  of  free-air  anomalies  in  pronounced  mountainous  terrain 
can  only  be  achieved  by  a  polynomial  of  high  degree  by  means  of  Ag-data 
available  in  a  network  of  high  resolution.  As  a_  consequence,  satisfac¬ 
tory  linear  interpolation  requires  small  mesh  sizes  Ax,  Ay.  It  is 
possible  to  write  the  gravity  anomaly  as 

Ag  »  Ag  +  C  *  r  '  (40) 

where  Ag^  is  the  isostatic  anomaly  valid  for  the  compensated  geoid, 
represents  the  aggregate  of  terms  computable  from  the  known  topography, 
and  r  is  a  random-type  error.  In  a  more  general  form,  also  applicable 
to  the  optimal  estimation  of  vertical  deflections,  eq.  (40)  is 
reformulated  as 

m  ■  s  +•  n  r  •  '  .  (41) 


3*  For  details  see  Baussus  von  Luetzow  [io] 
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In  this  equation,  m  is  a  "message"  variable,  s  is  a  "signal"  variable,  n 
is  deterministic  or  computable  "noise,"  and  r  is  random-type  noise. 

Under  consideration  of  a  linear  signal  estimation  structure,  a 
signal  can  then  be  optimally  estimated  as 

sQ  -  LCa^  -  n±  -  r±)  (42) 

where  L  denotes  a  linear  operator  and  the  subscripts  e  and  i  refer  to 
the  estimation  point  ?e  and  measurement  points  P^,  respectively.  The 
optimal  measurement  at  P9  results  in 

A  A 

m  *s  +  n  +r  =  L(m.  -  n.)  +  r  -  L  (r.)  +  n  .  (43) 

eeee  xx  e  1  e  '  " 

The  estimation  error  is 

e  a 

e(m  )  -  e(s  )  +  e  [r  -  L(r  )]  (44) 

9  9  9  1 

The  corresponding  estimation  error  resulting  from  the  utilization  of 
topographically  unmodified  measurements  m^  is 

e(m  )  *  s  -  L(s  )  +  r  -  L(r  )  +  n  -  L(n  ) 
e  e  x  e  1  e  1 

*  e(s@)  +  e[r@  -  L(ri)]  +  n@  -  L(ni).  (45) 

Comparison  of  eq.  (45)  with  eq.  (44)  shows  that  the  non-optimal  interpo¬ 
lation  process  is  associated  with  a  "topographic”  estimation  error  ne  - 
L(n^)  which  becomes  in  general  intolerable  in  moderate  to -rough 
mountainous  terrain  and  thus  induces  the  requirement  of  a  fine  mesh  data 
grid. 


Baussus  von  Luetzow  MoJ  obtained  the  following  solution  for  verti¬ 
cal  deflections: 

(CL  R  ff  \  .  “  „  ,A\  Jobs  a).  dS(»)  .  &g  +  SlUg)  /tn  Bn 

4mG  ^  8  Jl  ®  -^sin  a f  di|»  0  5  Itn  82^ 


II  LUi  *««  *  *&} 


C-6g2  -  Gi(6g3).  /tn  81 


/tn  sn 
I  tn  82  ' 


a 

In  this  equation,  0  denotes  global  mean  gravity,  &g  is  a  signal  variable 
profitably  synonymous  with  Ag^  in  eq.  (40),  a  is  the  azimuth  angle 
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counted  clockwise  from  north,  ij>  is  the  angle  between  the  radius  vectors 
r^  and  Tp  originating  at  the  earth's  spherical  center  in  accordance  with 
figure  2,  S  ($)  is  Stokes'  function,  and  62  are  northern  and  eastern 
terrain  inclinations,  respectively,  G^  is  a  specified  integral  function, 
5g1 ,  6g2  dg*  are  quantities  computable  from  the  earth's  topography,  C 
is  the  terrain  correction,  and  a  indicates  the  unit  sphere. 


Figure  2 


The  first  two  terms  of  eq.  (16),  involving  the  anomaly  hg,  represent  the 
"signal"  components  of  C  and  n.  The  following  three  terms  constitute 
computable  topographic  "noise”.  Permitting  random- type  errors  r^  and 
,  eq.  (4.6)  can  be  written  in  analogy  with  eqs.  (40)  and  (41)  as 


The  numerical  determination  of  the  three 'topographi :  terras  of  eq.  (46) 
13  a  complex  task,  which  can,  however,  be  accomplished  without  inherent 
difficulties  by  means  of  high-speed  computers.  The  estimation  of  t  and 
p  •'ropi  a  3et  of  given  vertical  deflections  can  be  achieved  by  means  of 
spatial  covariance  functions.  It  is  clear  from  the  above  analysis  that 
the  optimal  estimation  method  outlined  £n  section  4  requires  the 
simultaneous  consideration  of  eq.  (47).  Otherwise,  the  linear  aggre- 
-'pvolving  Z  and  n  in  eqs.  (35)  and  (36)  are  to  be  estimated  with 
er*o>*  4‘-om  a  limited  set  of  £  and  n  in  order  to  obtain  solutions 
without  linear  aggregates.  Repetitive  surveys  would  reduce  the  effects 
0 f  instrument  errors.  In  summary,  the  determination  of  deflections  of 
■ne  vertical  in  strongly  mountainous  terrain  requires  a  much  increased 
computational  effort  to  obtain  minimum  errgr  accuracies.  Without  such 
effort,  there  is  a  degradation  in  accuracy.  __  t 

6,  UTILIZATION  OF  GENERATED  GRAVITY  ’/ECTOR  COMPONENT  INFORMATION. 

Rapid  economical  gravity  vector  component  determination  is  necessary, 
useful,  or  potentially  useful  for  the  following  purposes: 
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(a)  Transformation  of  astronomic  coordinates  and  astronomic  azi¬ 
muths  into  ideal  terrestrial-geodetic  coordinates  and  geodetic  azimuths, 
requiring  knowledge  of  absolute  deflections  of  the  vertical,  indispens- 
ible  for  long  range  ballistic  missiles. 

(b)  Establishment  of  a  three-dimensional  anomalous  gravity  poten¬ 
tial  T  or  analytical  upward  continuation  of  first  order  derivatives  of 
Ag,  n,  and  £  from  discrete  measurements  thereof  at  the  earth's  surface, 
to  extend  over  relatively  large  areas,  indispensible  for  long  range 
ballistic  missiles  with  inertial  guidance. 

(c)  Adjustment  of  geodetic  networks. 

(d)  Geophysical  prospecting  including  estimation  of  the  deriva¬ 
tive  a  201?”^  +  GH-1?  tn  4  -  G  (~  +  )  in  semi-flat  terrain 

or  o  x  o  y 

where  N  stands  for  geoidal  undulation  or  height  anomaly.^ 

(e)  Determination  of  geoidal  differences  on  land. 

(f)  Construction  of  regional  gravity  anomaly  and  vertical  deflec¬ 
tion  networks  with  mesh  sizes  Ax  *  Ay  of  about  5  Km  as  input  information 
for  highly  accurate  gravity  programmed  inertial  surveying  systems. 

(g)  Accuracy  improvement  of  medium  -ange  ballistic  missiles  in 
conjunction  with  meteorological  corrections  and  of  cruise  missiles  by 
utilization  of  vertical .deflection  data  in  accessible  areas. 

7.  CONCLUSIONS.  The  determination  of  vertical  gravity  vector  compo¬ 
nents  or  gravity  anomalies  by  means  of  inertial  and  boundary  data  has 
been  very  promising  as  to  rapidity  and  accuracy.  High  accuracy  estima¬ 
tion  of  horizontal  gravity  vector  components  or  deflections  of  the 
vertical  from  inertial  and  boundary  data  requires  the  installation  of 
gryoscopes,  accelerometers,  and  velocity  quantizers  with  small  error 
variances  and  short  correlation  times,  the  development  of  sophisticated 
estimation  methods  consistent  with  an  outlined  structure,  and,  probably, 
a  three-dimensional  Kalman  filter.  Deflection  accuracy  in  strongly 
mountainous  terrain  is  degraded.  For  this  reason,  a  maximum  rms  error 
requirement  of  about  0.3  arcsec  in  both  semi-flat  and  mountainous  ter¬ 
rain  may  only  be  achieved  by  a  gradioraeter-aided  inertial  system.  Such 


4.  See  Telford  et  al  [ll],  pp.  7-104,  for  geophysical  applications  and 
Heiskanen  and  Moritz  [12],  p.  117,  for  the  equation. 
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a  system  would  al30  be  beneficial  for  the  location  of  smaller  subterra¬ 
neous  mass  anomalies.  An  effective  estimation  of  gravity  anomalies  and 
deflections  of  the  vertical  in  mountainous  terrain  from  data  at  points 
separated  by  distances  of  the  order  5-10  Km  is  possible  with  consid¬ 
erable  computational  effort  under  utilization  of  the  developed  equations 
and  employment  of  spatial  covariance  functions.  Numerous  military  and 
non-military  applications  of  gravity  vector  component  information  gener¬ 
ated  by  inertial  and  auxiliary  data  do  exist.  Gradiometer-aided  iner¬ 
tial  systems  offer  additional  advantages. 
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